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The van der Waals (VDW) equation of state predicts the existence of a first-order liquid-gas phase 
transition and contains a critical point. The VDW equation with Fermi statistics is applied to a 
description of the nuclear matter. The nucleon number fluctuations near the critical point of nuclear 
matter are studied. The scaled variance, skewness, and kurtosis diverge at the critical point. It is 
found that the crossover region of the phase diagram is characterized by the large values of the scaled 
variance, the almost zero skewness, and the significantly negative kurtosis. The rich structures of 
the skewness and kurtosis are observed in the phase diagram in the wide region around the critical 
point, namely, they both may attain large positive or negative values. 
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I. INTRODUCTION 

The first-order liquid-gas phase transition is a well 
known phenomenon that takes place in atomic and/or 
molecular systems, in the system of interacting nu¬ 
cleons (nuclear matter), and most probably between 
hadrons and quark-gluon plasma at large baryonic 
densities. In all these cases the phase transition line in 
the plane of temperature T and chemical potential p 
has an end point. It is called the critical point (CP), 
and it demonstrates some universal features typical 
for the second-order phase transitions, in particular, 
anomalously large fluctuations. 

The particle number fluctuations are characterized 
by the central moments, ((AN) 2 }, ((AN) 3 }, and 
((AN) 4 }, etc., where (...) denotes the event-by-event 
averaging and AN = TV — (TV). The scaled variance 
w[TV], skewness Sa, and kurtosis kct 2 defined as the 
following combinations of the central moments, 

((AAQ 2 ) 

(TV) ’ 

((ATV) 3 ) 

((AN) 2 } ’ 1 ’ 

((ATV) 4 ) - 3((ATV) 2 ) 2 
((ATV) 2 ) 

are the well known size-independent (intensive) mea¬ 
sures of particle number fluctuations. 

In the grand canonical ensemble (GCE) the pres¬ 
sure p plays the role of the thermodynamical potential, 
and its natural variables are temperature T and chem¬ 
ical potential /i. The particle number fluctuations can 


w[TV] = 
Sa = 


be characterized by the following dimensionless cumu- 
lants (n = 1,2,...), 
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The fluctuation measures in Eq. © can be then pre¬ 
sented as the following: 




( 3 ) 


The study of event-by-event fluctuations in high- 
energy nucleus-nucleus collisions opens new possibil¬ 
ities to investigate properties of strongly interacting 
matter (see, e.g., Refs. flj and {§] and references 
therein), and the experimental search for the chiral 
CP is now in progress (see, e.g., Refs. M and refer¬ 
ences therein). The fluctuation signals of the QCD CP 
were discussed in Refs. @-@|. In particular, the non- 
Gaussian fluctuation measures of conserved charges 
such as the skewness Sa and kurtosis kct 2 have at¬ 
tracted much attention recently (see, e.g., Refs, [sj] 
and 0 ). The higher moments of conserved charges 
were suggested as probes to study the QCD phase 
structure 00 , and have been calculated in various 
effective QCD models 0-0- Experimentally, the 
higher moments of net-proton and net-charge multi¬ 
plicity were recently measured by the STAR collabo¬ 
ration in Au+Au collisions in ^/s NN = 7.7 — 200 GeV 
energy range 00. However, no definitive conclu¬ 
sion regarding the existence and location of the chiral 
CP has been obtained yet. 

In the present paper the scaled variance, skewness, 
and kurtosis of net-nucleon (net-baryon) number fluc¬ 
tuations near the critical point of nuclear matter are 
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studied. A presence of the liquid-gas phase transi¬ 
tion in nuclear matter was reported in a large number 
of papers, see, e.g., Refs. [2014221 ] . Experimental esti¬ 
mates of the nuclear matter CP, T c = 17.9 MeV and 
n c = 0.06 fm -3 , were presented recently in Ref. [0 . 
At such small temperatures the effects of deconfine¬ 
ment and of production of new particles, such as pi- 
ons, are expected to be negligible, and the number 
of nucleons is essentially a conserved quantity. Thus, 
very different physical pictures of the critical behav¬ 
ior in nuclear matter and near the chiral CP are evi¬ 
dently expected. Nevertheless, in both cases the fluc¬ 
tuations of conserved charges are expected to be sen¬ 
sitive probes of critical behavior, and may be used to 
pinpoint the location of the corresponding CP. 


II. NUCLEAR MATTER WITH THE VAN 
DER WAALS EQUATION OF STATE 

In this work we use the van der Waals (VDW) equa¬ 
tion of state to study the measures © of particle num¬ 
ber fluctuations near the CP of the nuclear matter. 
The VDW model contains the first-order liquid-gas 
phase transition which ends at the CP. In the canoni¬ 
cal ensemble the classical VDW equation of state has 
a simple and transparent form (see, e.g., Ref. 0 ): 


the system entropy. The generalization of the VDW 
equation which includes effects of the quantum statis¬ 
tics was recently proposed in Ref. [2fij . The pressure 
and the particle density are then defined by the fol¬ 
lowing system of two equations for p(T, p) and n{T , p) 
functions: 


where 


p(T,p) =p ld (T,p*)-a n 2 (T, p ), 
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p* = /i — bp(T, p) — abn 2 (T, p) + 2 an(T, p) . (8) 

The p ld and n ld are the expressions for the quantum 
ideal gas pressure and particle density, respectively: 
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where V is the system volume, a > 0 and b > 0 are 
the VDW parameters that describe attractive and re¬ 
pulsive interactions, respectively, and n = N/V is the 
particle number density. The CP corresponds to the 
temperature T c and particle number density n c , where 


dp 

dn 


T 


= 0, 



(5) 


In order to apply the VDW equation of state to 
systems with a variable number of particles the GCE 
formulation is needed (see Ref. 0). 

In the following we consider the VDW equation of 
state for the system of interacting nucleons. We re¬ 
strict our consideration to temperatures T < 40 MeV, 
thus, the production of new particles (like pions) is 
neglected. In addition, both the nucleon clusters (i.e., 
ordinary nuclei) and the baryonic resonances (like N* 
and A) are neglected. Within these approximations, 
the number of nucleons, N , becomes a conserved num¬ 
ber, and the chemical potential p of the GCE regulates 
the number density of nucleons. 

At low temperatures and/or high particle number 
densities the Boltzmann approximation becomes in¬ 
adequate and leads to unphysical negative values of 


where g is the degeneracy factor (the number of spin 
and isospin states) and m is the particle mass. In 
Eqs. © and (flOl) . rj = +1 for Fermi statistics, rj = 
— 1 for Bose statistics, and p = 0 for the Boltzmann 
approximation. 

The VDW equation of state with Fermi statistics 
was used in Ref. 0 to describe the properties of sym¬ 
metric nuclear matter. In this case, p = 1, g = 4, and 
m = 938 MeV in Eqs. © and (flTill . Parameters a 
and b are fixed to reproduce the properties of nuclear 
matter in its ground state, i.e., at T = 0 it should be 
p = 0, n = no — 0.16 fnU 3 , and the binding energy 
per nucleon Eb = E/N — m = —16 MeV. These con¬ 
ditions give a = 329 MeV fm 3 and b = 3.42 fm 3 . Note 
that particle volume parameter b is connected to its 
hard-core radius r = [36/(167T)] 1 / 3 = 0.59 fm. In the 
GCE, at fixed T and p, equations © and © may 
have more than one solution. In such a case a solu¬ 
tion with the largest pressure is selected in accordance 
with the Gibbs criteria (see Ref. [0 for details). 


III. NUCLEON NUMBER FLUCTUATIONS 
NEAR THE CRITICAL POINT 

The phase transition line, p = p m i x (T), shown in 
Fig. m (a), starts from the normal nuclear matter 
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Figure 1: (Color online) The (a) particle number density n(T, /x) © and (b) scaled variance ux[7V] (HID calculated for 
the symmetric nuclear matter in ( T , /x) coordinates within VDW equation of state for fermions. The open circle at T = 0 
denotes the ground state of nuclear matter, the solid circle at T = T c corresponds to the CP, and the phase transition 
curve fi = /Xmix(T) is depicted by the solid line. 


state with T = 0, fio == 922 MeV and ends at the 
CP with T c = 19.7 MeV, /x c = 908 MeV (this gives 
n c = 0.07 fm -3 ) 1 . At each point of the phase transi¬ 
tion line, two solutions with different particle densities 
(the liquid and gas states) and equal pressures exist, 
i.e., this is a line of the first-order phase transition. 
At T > T c only a single solution n(T, /x) exists. Nev¬ 
ertheless, as seen from Fig. |T] (a), a rapid although 
continuous change of particle number density takes 
place in a narrow T-fi region (the so-called crossover 
region) even at T > T c . 

Using Eq. © one calculates the scaled variance 
cu[V] as 


k 2 T (dn 
" [JV| = h = n U; 
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1 The Boltzmann approximation rj = 0 leads to n c = 1/3 b = 
0.10 fm” 3 and T c = 8a/(27b) “ 28.5 MeV. 


with 77 = 1 corresponds to the scaled variance of par¬ 
ticle number fluctuations in the ideal Fermi gas (in 
the Boltzmann approximation, 77 = 0, it is reduced to 
td id — 1)- 

It is clearly seen from Eq. (fTTl) that the repulsive in¬ 
teractions suppress the particle number fluctuations, 
whereas the attractive interactions lead to their en¬ 
hancement. The scaled variance dill) is shown in Fig. [T] 
(b). At any fixed value of temperature, uj[N] —> 1 as 
/i decreases. In this case, n —> 0 and uiid(T, fx*) —> 1, 
thus, the Boltzmann ideal gas results are recovered. 
The scaled variance becomes small, w[fV] -C 1, as fi 
increases. In this case, the particle number density 
goes to its limiting value, n —> 1/b. The scaled vari¬ 
ance m diverges at the CP (note that the thermody¬ 
namic limit V —y 00 is assumed). As seen from Fig. [T] 
(b) the large values of w[V] 1 take place along the 
crossover region, even far away from the CP. 

Using Eq. © one also calculates the skewness 




du[N] 

dfx 


and the kurtosis 
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k 2 
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(13) 


(14) 


shown in Figs. [ 2 ] (a) and (b), respectively. Similarly to 
the scaled variance, the skewness and kurtosis diverge 
at the CP. However, these higher moments of the par¬ 
ticle number distribution show much richer structures: 
the behavior of Sa and na 2 crucially depend on the 
path of approach to the CP. 
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Figure 2: (Color online) The (a) skewness Sa (1131) and (b) kurtosis na 2 (11-11) calculated for the symmetric nuclear 
matter in ( T , /x) coordinates within VDW equation of state for fermions. 


As seen from Fig. [2] (a) the liquid phase corresponds 
to Sa < 0, whereas the gas phase corresponds to 
Sa > 0. The line of S'cr = 0 goes from the critical 
point along the crossover region, and the T-/x regions 
with Sa 1 and Sa -C —1 are placed just under and 
above this line, respectively. From Fig. [2] (b) it is also 
seen that the crossover region of the phase diagram 
is characterized by the significantly negative kurtosis, 
na 1 2 3 4 < -1. However, outside this crossover region one 
observes large positive values of the kurtosis, na 2 1, 
in a rather wide T-/x area around the CP of nuclear 
matter. Qualitatively, our findings with regards to the 
fluctuation patterns near the CP are consistent with 
previous results based on effective QCD models (see, 
e.g., Refs, mmmi or on model-independent uni¬ 
versality arguments with regards to critical behavior 
in the vicinity of the QCD critical point 

IV. SUMMARY 

In summary, the fluctuation signatures of the nu¬ 
clear matter CP — scaled variance, skewness, and 
kurtosis — are calculated within the quantum formu¬ 


lation of the VDW equation of state. The scaled vari¬ 
ance diverges, a;[AT] —► oo, at the CP, and the large 
values of w[iV] 1 take place along the crossover re¬ 
gion even far away from the CP. The behavior of Sa 
and na 2 at the CP is more complicated. The limit¬ 
ing singular values of these quantities depend on the 
path of approach to the CP. The rich structures of 
the skewness and kurtosis are observed in the wide 
T-/i area around the CP of nuclear matter. We hope 
that results obtained in this paper can be useful for 
the identification of the CP signatures of nuclear mat¬ 
ter in heavy-ion collision experiments. 
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